Effect of the Seasonal Climatic Variations on the Flavonoid Accumulation in Vitis vinifera cvs. ‘Muscat Hamburg’ and ‘Victoria’ Grapes under the Double Cropping System

Under the double cropping system, berries usually showed significant quality variations in the summer and winter seasons. In the two-year/four-consecutive-season study, two table grapes of ‘Muscat Hamburg’ and ‘Victoria’ were investigated to determine the phenolic compounds in their berries. Different from those of the summer season, the berries in the winter season suffered no high-temperature stress since veraison to harvest in 2014 and 2015. The variations in the season temperatures led to a higher anthocyanin concentration in the winter season berries of ‘Muscat Hamburg’ grapes than that in the summer berries, while the summer season berries had higher proportions of acylated and methylated anthocyanins than those in the winter season berries. Similar to the anthocyanins, the winter season berries also had a higher flavonol concentration in both varieties. Transcriptome analysis showed that the upregulated genes involved in the flavonoid pathway in the winter season berries were agreed with the changes found in the metabolites. However, the influence of the growing seasons on the flavanols was not consistent in the two varieties, and the variations in VviLARs between the grapes of ‘Muscat Hamburg’ and ‘Victoria’ might be the cause. This research helped us better understand the double cropping system and how the climate factors affected the phenolic compounds in the double cropping system.


Introduction
In South China, the subtropical humid monsoon climate was usually considered suboptimal for grape cultivation because of the extremely high temperature and the concentrated rainfall in the growing season, especially in summer [1]. However, table grapes had a rapid development in South China in recent years, such as in Guangxi province. Cheng et al. [2] showed that the viticultural area of Guangxi had increased by three times, and the annual production value increased from 246 million to 2.6 billion yuan (CNY) since the beginning of this century. The rain shelter cultivation and double cropping system techniques contributed to the booming development. The double cropping system was a common technique in the subtropical regions [3,4]. In these regions, the excessive rainfall in the summer led to fungal diseases and increased the rot incidence, which negatively affected the grape quality [4,5]. After applying the double cropping system, berries ripened earlier than those with the traditional single cropping system, making it possible to avoid the intense rainfall and heatwave in the summer season. Furthermore, berries in the winter season usually had a better quality than those in the summer season [6,7].
Flavonoid compounds are critical secondary metabolites in grapes, including anthocyanins, flavonols, flavanols, etc. [8]. These compounds are synthesized through the phenylpropanoid-flavonoid pathway and share the common upstream steps through phenylalanine ammonia-lyase (PAL) to flavanone 3-hydroxylase (F3H) [1]. The accumulation of flavonoids is sensitive to environmental changes such as temperature and light [9]. Previous studies showed that the relatively low temperature was beneficial for the accumulation of anthocyanins, while high temperature inhibited the synthesis of anthocyanins [10,11]. Cohen et al. [12] showed that diurnal temperature variation altered the initial rates of proanthocyanidin accumulation, which was correlated strongly with the expression of the core genes in the flavonoid pathway. As for the light effect, the increased cluster exposure could promote the synthesis of anthocyanins and flavonols, while shading treatment would inhibit the related metabolism of these compounds [13,14]. In addition, changes in the climate factors could also affect the portion of various flavonoids. For example, the high temperature would increase the proportion of coumalylated anthocyanins in the 'Merlot' grape [15]. Water deficit could increase the proportion of the 3 5 -hydroxylated and the methoxylated anthocyanins in the 'Cabernet Sauvignon' grape [16].
In the present study, we chose two table grapes that occupy a more extensive market than the wine grapes in South China to determine the variations in the flavonoid profiles. The 'Victoria' is a table grape variety cultivated in Puglia region at the Horticulture Research Institute of Bucharest [17]. The 'Muscat Hamburg' is a classical cultivar of black table grape grown in many parts of Europe, highly appreciated for its beautiful bunches and pleasant Muscat flavor [18]. Although many studies focus on the aroma characteristics of the two varieties, their phenolic profiles under the double cropping were few reported. There were significant climate variations between the summer and winter seasons under the double cropping system, which usually exceeded the effects of the vintages in the traditional viticultural regions. The results could help the researchers better understand how climate parameters affected grape phenolic compounds under the double cropping system, which might improve the understanding of the grape berries in response to the climate changes accompanied by extreme weather conditions in the future

Vineyard and Double Cropping System
Field experiments were performed in V. vinifera L. cvs. 'Muscat Hamburg' (MH) and 'Victoria' (V) grapevines (both grafted on 'SO4 ) in 2014 and 2015 in four consecutive growing seasons. Experiment site located in Guangxi Academy of Agricultural Sciences which was in South China (22 • 36 N-108 • 14 E, elevation 104 m). Vines were planted in 2007 and trained to a Y-shaped training system with 2 × 4/5 shoots per meter and 1.0 m cordon above ground. Rain-shelters were applied to all vines to prevent over-rainfall damage. The vineyard was north-south row orientation with inter-and intra-row spacing of 3.5 m × 1 m. The double cropping system in the experiment site was described by Chen et al. [1]. In each vintage, vines had two growing seasons: the summer season and the winter season. The summer season began in March and the grapes were harvested around late July and early August. Then vines were pruned and enforced with 2.5-3.0% hydrogen cyanamide to start the winter season. Vines were budburst in August and grapes were harvested in January next year.

Berry Sampling and Meteorological Data Collection
In each growing season, berries were sampled at four E-L stages [19]: (1) pea-size (E-L 31), (2) onset of veraison (E-L 35), (3) veraison complement (E-L 36) and (4) harvest (E- L 38). There were three biological replicates for each variety. For each replicate, 300 berries were randomly sampled from about 50 vines. Berries were immediately frozen in liquid nitrogen and stored at −80 • C for metabolite and transcriptome analysis. Meteorological data was acquired from a local climate monitoring station within one kilometer away from the experiment site. Photosynthetically active radiation and temperature were recorded per hour. Accumulated rainfall was recorded per day.

Extraction of Grapes Phenolic Compounds
For each replicate, 100 berries were selected randomly to peel off the skins in frozen status. Skins were quickly frozen in liquid nitrogen and then grounded to powder. The powder was dried at −40 • C under vacuum conditions for 24 h.
The extraction of anthocyanins and flavonols in skins was according to a previous study [20]. Briefly, 0.1000 g (±0.0002 g) dried powder was accurately weighed and put into the centrifuge tube. Then the dried powder was macerated and sonicated in 50% (v/v) methanol in water (1.0 mL) for 20 min under a low temperature (4 • C) and dark condition. The mixture was centrifugated for 10 min at 8000 rpm to acquire the supernatant. The residues were re-extracted again and all the supernatants were combined.
The extraction of flavanols in skins was according to Liang et al. [21]. Briefly, 0.1000 g (±0.0002 g) dried powder was accurately weighed and put into the centrifuge tube. For the various flavan-3-ol units, 1 mL of phloroglucinol buffer (0.5% ascorbate, 300 mmol/L HCl and 50 g/L phloroglucinol in methanol) was used to mix the dried powder and incubated at 50 • C for 20 min. Then the mixture was neutralized with 1.0 mL aqueous sodium acetate (200 mM) and centrifuged for 15 min at 8000 rpm. The residue was extracted twice and all the supernatants were collected and combined. For the free flavan-3-ol monomers, 1 mL of 70% acetone with 0.5% ascorbate was used to mix the dried powder and then centrifuged for 15 min at 8000 rpm. The extraction of the residue was conducted twice and all the supernatants were collected and combined. Then the supernatants were dried using a nitrogen stream at 30 • C. The dried samples were dissolved in 200 µL acidified methanol with 1% (v/v) HCl and then neutralized with 200 µL aqueous sodium acetate (200 mM).

HPLC-MS Analysis of Phenolic Compounds in Grapes
Anthocyanins were analyzed with an Agilent 1100 series HPLC-MSD trap VL equipped with a diode array detector (DAD) and a Kromasil C 18 column (250 × 4.6 mm, 5 µm). Mobile phase A was 2% formic acid in water and B was 2% formic acid in acetonitrile. The detailed LC procedures and MS conditions have been described by He et al. [22]. Anthocyanins were quantified using the malvidin-3-O-glucoside as the external standards and expressed as mg/kg berry fresh weight (FW). Flavonols were analyzed on an Agilent 1200 series HPLC-MSD trap VL equipped with a variable wavelength detector (VWD) and a Zorbax EclipseXDB-C 18 column (250 × 4.6 mm, 5 µm). Mobile phase A consisted of acetonitrile/formic acid/water (50:85:865, v/v/v), and mobile phase B consisted of acetonitrile/methanol/formic acid/water (250:450:85:215, v/v/v/v). The detailed LC procedures and MS conditions were described by Chen et al. [1]. Flavonols were quantified using the quercetin-3-O-glucoside as the external standards and expressed as mg/kg berry fresh weight (FW). Flavanols were analyzed on an Agilent 1200 series HPLC system equipped with a Poroshell 120 EC-C18 column (150 × 2.1 mm, 2.7 µm) and Agilent 6410 QqQ instrument equipped with an electrospray ionization source. Mobile phase A was aqueous 0.1% formic acid and B was a mixture of acetonitrile/methanol (50:50, v/v) containing 0.1% formic acid. The detailed LC procedures and MS conditions have been described by Chen et al. [1]. Flavanols were quantified by using (+)-catechin, (-)-epicatechin, (-)epicatechin-3-O-gallate and (-)-epigallocatechin as the external standards and expressed as mg/kg berry fresh weight (FW).

RNA Extraction and Transcriptome Sequencing
Berries from three sampling points were selected (E-L 35, 36, and 38) for the RNA extraction. In this case, 50 berries were randomly selected from each biological replicate. Then the berries were de-seeded and smashed into powder under liquid nitrogen protection. The SpectrumTM Plant Total RNA Kit (Sigma-Aldrich, Carlsbad, CA, USA) was used for  Transcriptome analysis was conducted on the Illumina HiSeqTM  2000 platform with 50-bp single reads and aligned against the reference grapevine genome 12 × V2, allowing no more than two mismatches. The longest transcript was chosen to calculate the fragments per kilobases per million reads (FPKM) value when more than one transcript was obtained for a single gene. The R package 'DESeq2 was used to identify differentially expressed genes (DEGs), and the criteria were set as false discovery rate ≤ 0.05 and fold change ≥ 2. The data have been deposited in the NCBI Gene Expression Omnibus (GEO) database and are accessible through GEO accession GSE168785.

Statistical Analysis
The SPSS version 22.0 was used for all significance analysis at p < 0.05 (Duncan's multiple range test or t-test). The figures were prepared by using GraphPad Prism 8.0.2 (GraphPad Software, San Diego, CA, USA) and R statistical environment (3.6.1). Heatmap was prepared using the 'pheatmap' package in R. Principal component analysis (PCA) and orthogonal partial least-squares discrimination analysis (OPLS-DA) were performed in SIMCA 14.1 (Umetrics, Umea, Sweden).

Meteorological Data
The climate conditions of each development stage in 2014 and 2015 were shown in Table 1. Stage I was from the full bloom to the veraison beginning. Stage II was from the veraison beginning to the veraison completement. Stage III was from the veraison completement to the harvest. The 'Muscat Hamburg' and 'Victoria' grapes had similar phenological stages in the winter seasons of 2014 and 2015 (Supplementary Table S1). So they experienced similar climate conditions in these seasons. In the summer season of 2015, the 'Muscat Hamburg' grapes were harvested 23 days later than the 'Victoria' grapes. In the summer season, the GDD was significantly higher than in the winter season. The average daily temperature in the summer season was about 30 • C in three development stages. While in the winter season, only stage I could be up to 23 • C in terms of average daily temperature. In stage III of the winter season, the average temperature was no more than 16 • C, which indicated a cool weather station during the ripening period. Notably, the high-temperature hours in the summer season were at least three folds higher than those in the winter season. Furthermore, during stage II and stage III, almost no high-temperature weather occurred in the winter season. For the cumulative PAR/sunshine hours and the rainfall, there were no consistent trends in 2014 and 2015. In 2014, the winter season had more cumulative PAR/sunshine hours and rainfall than the summer season. However, in 2015, the winter season had fewer cumulative PAR/sunshine hours and rainfall than the summer season.
The proportion of the acylated and the B-ring substituted anthocyanins was calculated and shown in Figure 1. For acylated anthocyanins, the summer season berries had higher proportions of the acetylated and the coumarylated anthocyanins than the winter season berries in both of the two vintages. The proportion of the acetylated and the coumarylated anthocyanins increased consistently during the growing season in the summer season berries while keeping stable in the winter season berries. As for the methylated anthocyanins, the summer season berries also had a higher proportion than the winter season berries. The increased proportion of malvidin-based anthocyanins and decreased proportion of cyanidin-based anthocyanins caused a higher proportion of the methylated anthocyanins in the summer season berries. In terms of the 3 5 -hydroxylated anthocyanins, their proportions followed an increasing trend in all the growing seasons. Similar to acylated and methylated anthocyanins, the proportion of 3 5 -hydroxylated anthocyanins was also higher in the summer season berries than in the winter season berries. The 3 5 -hydroxylated anthocyanins consist of the malvidin-based, the delphinidinbased, and the petunidin-based anthocyanins. Among the three types of anthocyanins, the increased proportion of the malvidin-based anthocyanins caused a higher proportion of the 3 5 -hydroxylated in the summer season berries. season berries. The increased proportion of malvidin-based proportion of cyanidin-based anthocyanins caused a highe anthocyanins in the summer season berries. In terms of th nins, their proportions followed an increasing trend in all t acylated and methylated anthocyanins, the proportion of 3 was also higher in the summer season berries than in the w hydroxylated anthocyanins consist of the malvidin-based, t petunidin-based anthocyanins. Among the three types of an portion of the malvidin-based anthocyanins caused a hig droxylated in the summer season berries.

Flavonol Composition
As shown in Table 3

Flavonol Composition
As shown in Table 3   concentration among all the four seasons in the two varieties.

Flavanol Composition
Respect to flavonols in the grape skins, the basic monomer (+)-catechin, (-)-epicatechin, (-)-epicatechin-3-O-gallate and (-)-epigallocatechin were detected by the HPLC-MS in three forms: terminal subunits, extension subunits and free monomers (Table 4). Among the four basic monomers, (-)-epicatechin occupied the highest proportion (approximately 70%) of the total flavanol concentration. There were no consistent trends in the (-)-epicatechin concentration between the summer and winter season berries in the two varieties. In 2014, the winter season berries had a higher (-)-epicatechin concentration than the summer season berries in the 'Muscat Hamburg' grape, while had an opposite result in the 'Victoria' grapes. However, in 2015, there was no significant difference in (  two varieties. In 2014, the winter season berries had a higher (-)-epicatec than the summer season berries in the 'Muscat Hamburg' grape, while result in the 'Victoria' grapes. However, in 2015, there was no significant epicatechin concentration between the summer and the winter season compounds, only (-)-epicatechin-3-O-gallate showed significant differe summer and the winter season berries in the two varieties. The summ had a higher (-)-epicatechin-3-O-gallate concentration than the winter 2014 but the opposite result showed in 2015. In terms of total flavanol the summer and the winter season berries (Figure 3), they showed decre the development stages which peaked at E-L 31. In the 'Muscat Hamburg winter season berries had the highest flavanol concentration among all at harvest. While in the 'Victoria' grapes, the 2014 summer season berrie flavanol concentration among all the seasons at harvest. Values are reported as means ± SD of three biological replicates. Different letter indicate significant differences among treatments (Duncan's multiple range test

Principal Component Analysis (PCA) and Orthogonal Partial Least-Squa Analysis (OPLS-DA) Based on the Phenolic Profiles at Different Stages
To better understand how the phenolic profiles varied in 'Muscat H and 'Victoria' grapes, as well as in different seasons. The principal com (PCA) was used to classify the different samples, which consisted of all ment stages of two varieties in different seasons, as shown in Figure 4 an Figure S1. The first two principal components explained 70.4% of the to accounted for 51.8% of the total variance, which could separate samples

Principal Component Analysis (PCA) and Orthogonal Partial Least-Squares Discrimination Analysis (OPLS-DA) Based on the Phenolic Profiles at Different Stages
To better understand how the phenolic profiles varied in 'Muscat Hamburg' grapes and 'Victoria' grapes, as well as in different seasons. The principal component analysis (PCA) was used to classify the different samples, which consisted of all the four development stages of two varieties in different seasons, as shown in Figure 4 and Supplementary Figure S1. The first two principal components explained 70.4% of the total variance. PC1 accounted for 51.8% of the total variance, which could separate samples from MH and  (Figure 1a). In the loading plot (Figure 1b), PC1 was characterized by all the anthocyanins and some of the flavonols. The V grapes had no anthocyanins so all the anthocyanins were located on the positive of axis x, which was in agreement with the previous analysis. Furthermore, only quercetin and kaempferol-based flavonols were detected in the V grapes. Other types of flavonols were also located on the positive of axis x. Except for the variations in different varieties, the samples from different seasons and development stages were also marked in the PCA, as shown in Supplementary Figure S1. However, the PCA could not separate these samples clearly. So the variations in the variety characteristics beyond the effects of growing seasons and development stages on grape phenolics. However, with respect to the summer and winter season effects (Supplementary Figure S1b), the samples from summer seasons were close to the coordinate origin, while some samples from winter seasons were abundant in anthocyanins and flavonols.
Foods 2022, 10, x FOR PEER REVIEW cyanins and some of the flavonols. The V grapes had no anthocyanins so all the a anins were located on the positive of axis x, which was in agreement with the p analysis. Furthermore, only quercetin and kaempferol-based flavonols were det the V grapes. Other types of flavonols were also located on the positive of axis x for the variations in different varieties, the samples from different seasons and d ment stages were also marked in the PCA, as shown in Supplementary Figure S1 ever, the PCA could not separate these samples clearly. So the variations in the characteristics beyond the effects of growing seasons and development stages o phenolics. However, with respect to the summer and winter season effects (Supp tary Figure S1b), the samples from summer seasons were close to the coordinate while some samples from winter seasons were abundant in anthocyanins and flav To better discriminate the effect of different growing seasons on grapes pheno files. The OPLS-DA model was used which had a better focus toward the studied o than PCA, as shown in Figure 5. The model has passed 200 permutation tests ( mentary Figure S2), indicating good fit and predictive abilities. Results showed samples from the summer and winter seasons could be separated better than in mentary Figure S1b. In the S-plot (Figure 5b), most phenolic compounds were loc the first quadrant of the coordinate axis except for EGC. So it was confirmed that t ter berries had more abundant phenolic compounds irrelevant with varieties or ot tors. The peonidin-3-O-glucoside contributed the highest variation in all phenol pounds. In the winter berries, the peonidin-3-O-glucoside was at least 8 folds high in summer berries, which was shown in Table 2. For epigallocatechin, there was nificant difference between summer and winter berries. To better discriminate the effect of different growing seasons on grapes phenolic profiles. The OPLS-DA model was used which had a better focus toward the studied objective than PCA, as shown in Figure 5. The model has passed 200 permutation tests (Supplementary Figure S2), indicating good fit and predictive abilities. Results showed that the samples from the summer and winter seasons could be separated better than in Supplementary Figure S1b. In the S-plot (Figure 5b), most phenolic compounds were located on the first quadrant of the coordinate axis except for EGC. So it was confirmed that the winter berries had more abundant phenolic compounds irrelevant with varieties or other factors. The peonidin-3-O-glucoside contributed the highest variation in all phenolic compounds. In the winter berries, the peonidin-3-O-glucoside was at least 8 folds higher than in summer berries, which was shown in Table 2. For epigallocatechin, there was no significant difference between summer and winter berries.

Flavonoid Biosynthesis
To further understand the variation between the summer and the winter seas ries at the transcriptome level, the differentially expressed genes related to flavono synthesis were selected, as shown in Figure 6. The biosynthesis of flavonoids sha common upstream pathway through phenylalanine to dihydrokaempferol. So portant enzymes such as phenylalanine ammonia-lyase (PAL), cinnamate 4-hydr (C4H), 4-coumarate: CoA ligase (4CL), chalcone synthase (CHS), chalcone iso (CHI), flavanone 3-hydroxylase (F3H) were all involved in this part [24]. In the ' Hamburg' grapes, almost all the genes had higher expressions in the winter season than in the summer season berries in the upstream pathway. Even at E-L 38, thes still had higher expression in the winter season berries. In the 'Victoria' grapes, lected differentially expressed genes had higher expressions at E-L 35 and/or at while downregulated at E-L 38, such as VviPAL (VIT_216s0039g01100) and (VIT_211s0052g01090). Flavonoid 3′-hydroxylase (F3′H) and flavonoid 3′5′-hydr (F3′5′H) were involved in the two branch pathways of catalyzing the synthesis o substituted and the 3′5′-substituted flavonoids, respectively. There were two V (VIT_203s0063g01690 and VIT_209s0002g01090) selected as the differentially exp genes in the summer and the winter berries in the two varieties. Both of the two s VviF3′Hs had higher expression in the winter season berries, esp VIT_203s0063g01690, which had a higher expression at E-L 35, E-L 36 and E-L 38 F3′5′H, six VviF3′5′Hs were differentially expressed in the summer and the winter berries in the 'Muscat Hamburg' grapes, while the VviF3′5′Hs were almost not exp in the 'Victoria' grapes (Supplementary Table S2). Different from the red grape va the extremely low expressions of VviF3′5′Hs in the white grapes suggested that zyme flavonoid 3′,5′-hydroxylase was not expressed in their skins [23,25].

Flavonoid Biosynthesis
To further understand the variation between the summer and the winter season berries at the transcriptome level, the differentially expressed genes related to flavonoid biosynthesis were selected, as shown in Figure 6. The biosynthesis of flavonoids shared the common upstream pathway through phenylalanine to dihydrokaempferol. Some important enzymes such as phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), 4-coumarate: CoA ligase (4CL), chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H) were all involved in this part [24]. In the 'Muscat Hamburg' grapes, almost all the genes had higher expressions in the winter season berries than in the summer season berries in the upstream pathway. Even at E-L 38, these genes still had higher expression in the winter season berries. In the 'Victoria' grapes, the selected differentially expressed genes had higher expressions at E-L 35 and/or at E-L 36, while downregulated at E-L 38, such as VviPAL (VIT_216s0039g01100) and Vvi4CL (VIT_211s0052g01090). Flavonoid 3 -hydroxylase (F3 H) and flavonoid 3 5 -hydroxylase (F3 5 H) were involved in the two branch pathways of catalyzing the synthesis of the 3 -substituted and the 3 5 -substituted flavonoids, respectively. There were two VviF3 Hs (VIT_203s0063g01690 and VIT_209s0002g01090) selected as the differentially expressed genes in the summer and the winter berries in the two varieties. Both of the two selected VviF3 Hs had higher expression in the winter season berries, especially VIT_203s0063g01690, which had a higher expression at E-L 35, E-L 36 and E-L 38. As for F3 5 H, six VviF3 5 Hs were differentially expressed in the summer and the winter season berries in the 'Muscat Hamburg' grapes, while the VviF3 5 Hs were almost not expressed in the 'Victoria' grapes (Supplementary Table S2). Different from the red grape varieties, the extremely low expressions of VviF3 5 Hs in the white grapes suggested that the enzyme flavonoid 3 ,5hydroxylase was not expressed in their skins [23,25]. modification of anthocyanins have been widely studied, in (VIT_216s0039g02230), VviAOMT (VIT_201s0010g03510) (VIT_203s0017g00870) [28][29][30]. These genes were only expressed in ' grapes. The expressions of VviUFGT and VviAOMT were upregulated i berries, which was consistent with the coordinated expression of the u

Expression Profiles of Flavonoid Related Transcription Factors
In addition to the structural genes that encode enzymes in the f the related regulatory genes that control the transcription of these biosy also analyzed, as shown in Figure 7. Genes of the MYBA family were i ing the anthocyanin biosynthesis in the grapes via regulating the expr gene [31,32]. MYBA2 and MYBA3 were expressed in both 'Muscat Ha ria' grapes and upregulated in the winter season berries, while MY pressed in 'Muscat Hamburg' grapes and upregulated in the winter s Among the five known VviFLSs, only VviFLS4 (VIT_218s0001g03470) and VviFLS5 (VIT_218s0001g03430) were reported to express in grapes [26]. In the present study, VviFLS4 had a low expression in the 'Muscat Hamburg' grapes while it had a higher expression in the 2014 winter season of the 'Victoria' grapes (Supplementary Table S2). Compared to the summer season berries, the winter season berries had higher expressions of VviFLS4 and VviFLS5 than the summer season berries in both varieties. LAR and ANR were key enzymes in the production of the flavan-3-ol monomers required for the formation of proanthocyanidin polymers [27]. Two VviLARs and one VviANR were selected as the differentially expressed genes between the summer and the winter season berries. VviLAR1 (VIT_201s0011g02960) and VviLAR2 (VIT_217s0000g04150) were upregulated in the 'Muscat Hamburg' winter season berries. However, in the 'Victoria' grapes, VviLAR2 was downreg-ulated in the winter season berries. The genes involved in the synthesis and modification of anthocyanins have been widely studied, including VviUFGT (VIT_216s0039g02230), VviAOMT (VIT_201s0010g03510) and Vvi3AT (VIT_203s0017g00870) [28][29][30]. These genes were only expressed in 'Muscat' Hamburg grapes. The expressions of VviUFGT and VviAOMT were upregulated in the winter season berries, which was consistent with the coordinated expression of the upstream genes.

Expression Profiles of Flavonoid Related Transcription Factors
In addition to the structural genes that encode enzymes in the flavonoid pathway, the related regulatory genes that control the transcription of these biosynthetic genes were also analyzed, as shown in Figure 7. Genes of the MYBA family were involved in regulating the anthocyanin biosynthesis in the grapes via regulating the expression of the UFGT gene [31,32]. MYBA2 and MYBA3 were expressed in both 'Muscat Hamburg' and 'Victoria' grapes and upregulated in the winter season berries, while MYBA1 was only expressed in 'Muscat Hamburg' grapes and upregulated in the winter season berries. The regulators named VviMYBPA1, VviMYB5a, VviMYB5b and VviMYBC2-L1 were involved in the regulation of the proanthocyanidin-specific biosynthesis [33][34][35][36]. In the winter season berries, the expressions of VviMYB5a and VviMYB5b were upregulated in both varieties, while VviMYBPA1 was downregulated at E-L 36 in the 'Muscat Hamburg' grapes. Similarly, the VviMYBPA1 expression was also found to be upregulated with high temperatures [37], which was in agreement with our study. VviMYBC2-L1 was a negative repressor of the proanthocyanidin biosynthesis but showed opposite trends in the two varieties in terms of the winter vs. the summer season. The regulator of VviMYBF1 was involved in the regulation of flavonol biosynthesis, which could induce the expression of flavonol synthase (VviFLS1/VviFLS4), a key step of the initial flavonol pathway [38]. However, VviMYBF1 was downregulated in winter berries in both of 'Muscat Hamburg' and 'Victoria' grapes, which showed opposite trends to the upregulated VviFLS4 and higher flavonol concentration in the winter season berries. So there might be other transcription factors that play an important role in regulating VviFLS4. For example, VvibZIPC22 expression was induced by ultraviolet light (UV), accompanied by the expression of VviFLS4 and the accumulation of flavonols [39]. In the winter season berries, the expression of VvibZIPC22 was higher than that in the summer berries at E-L 35 and E-L 36, which might lead to a higher VviFLS4 expression.
VvibZIPC22 expression was induced by ultraviolet light (UV), acco sion of VviFLS4 and the accumulation of flavonols [39]. In the w expression of VvibZIPC22 was higher than that in the summer b 36, which might lead to a higher VviFLS4 expression.

Effects of Growing Season on Anthocyanin Composition of 'Musca
Under the double cropping system, the anthocyanin concent son berries was at least seven-fold higher than that in the summ study. The same result was also found by previous studies [6,7]. X the total anthocyanin content in 'Kyoho' skins in the winter season fold higher than the value of the corresponding summer season ported that anthocyanin content in the winter season berries of 'M be eleven-fold higher than that in the summer season season. The

Effects of Growing Season on Anthocyanin Composition of 'Muscat Hamburg' Grape
Under the double cropping system, the anthocyanin concentration in the winter season berries was at least seven-fold higher than that in the summer season berries in our study. The same result was also found by previous studies [6,7]. Xu et al. [7] showed that the total anthocyanin content in 'Kyoho' skins in the winter season berries was about fivefold higher than the value of the corresponding summer season berries. Zhu et al. [6] reported that anthocyanin content in the winter season berries of 'Muscat Hamburg' could be eleven-fold higher than that in the summer season season. The vast variation between the summer and the winter season berries in the anthocyanin accumulation showed a great range of phenotypic plasticity caused by climate factors. Among all climate factors, the high-temperature effect was the most conspicuous in our study. However, a different result was found by Chou and Li [3] that anthocyanin concentration in the 'Kyoho' grape was not affected by the seasonal variations between the summer and winter cropping system. They inferred that the complex environmental or physiological factors might overwhelm the influence of the temperature on the anthocyanin accumulation, although the temperature of the post veraison period was also higher in the summer cropping cycle than in the winter cropping cycle in their study [3]. Furthermore, the temperature variation between the seasons in their research was not notable, while in our study, the average daily temperature in the summer season could be 9 • C higher than that in the winter season, which should be resulted from the different climates between the studies. The expressions of VviUFGT and VviAOMT were upregulated in the winter season berries, which was consistent with the coordinated expression of the upstream genes in our study. In previous studies, the expression of the UFGT gene was significantly down-regulated by the high temperature until the mid-ripening in the 'Malbec' grapes [40] and from the mid-ripening to maturity in the 'Sangiovese' grapes [41]. In the winter season, almost no high-temperature weather occurred, which provided beneficial conditions for the anthocyanin biosynthesis.
Among all climate factors, the temperature should also be the dominant one in affecting the proportion of each group of anthocyanins. In a previous study, Tarara et al. [42] found that higher berry temperatures led to a higher proportion of the malvidin-based anthocyanins while decreasing the cyanidin-based anthocyanins, which was in agreement with our study. Acylated anthocyanins were known to be more stable than their nonacylated counterparts [43]. In the summer season, berries suffered more high-temperature stress than those in the winter, leading to the degradation of various anthocyanins types. Among all anthocyanins, non-acylated and non-methylated anthocyanins were easily degraded, leading to a higher proportion of the acylated and the methylated anthocyanins. The same result was also found by Tarara et al. [42]. They showed that grape berries might shunt more of the available anthocyanin toward acylation in response to the temperature stress, with the potential advantage to the plant of color stability within the vacuole because of the stability of the acylated compounds. As for the increased proportion of methylated and 3 5 -hydroxylated anthocyanins in the summer berries, the increased proportion of malvidin-based anthocyanins was the cause.

Effects of Growing Season on Berries Flavonol Composition
Flavonol was well known to be positively correlated with sunlight exposure. The biosynthesis of flavonol was upregulated by the solar radiation, leaving a fingerprint on the flavonol profile [44]. However, in our study, although the 2015 summer season had the most abundant sunshine hours and cumulative PAR, the flavonol concentration of the summer season berries was lower than that of the winter season berries in the 'Muscat Hamburg' grapes. The same results were also found by Chen et al. [1] and Zhu et al. [6] that the flavonol concentration in the winter season berries was higher than those in the summer season berries. Although flavonols responses to the temperature seemed to vary depending on the experimental parameters, the studies with whole-vine and the detached berry levels reported an effect on the flavonol concentration [41,45,46]. Pastore et al. [41] reported that in a low temperature condition, the flavonol concentration in berries was three times of that in the high-temperature conditions. In the summer season, the whole vine suffered more heat stress, thus causing a general decrease in metabolism at the whole vine level. The reduced primary metabolites as flavonoid precursors could impact the downstream production of flavonols [45]. So the decrease of the flavonol concentration in the summer season berries might mainly result from the high-temperature stress, which had a greater effect than the offset from the higher light radiation in the summer season of a certain vintage. However, compared to the summer (2014 summer vs. 2015 summer) or the winter (2014 winter vs. 2015 winter) growing seasons, more sunshine hours during the growing period were associated with the higher flavonol concentration in the berries in our study. Compared to the summer season berries, the winter season berries had higher expressions of VviFLS4 and VviFLS5 than the summer season berries in both of the two varieties. In general, the white grapes seemed to have a lower flavonol concentration than the red grapes. Mattivi et al. [23] showed that the total amount of flavonols found after the hydrolysis of the grape extracts ranged from 3.81 to 80.37 mg/kg, with a mean of 32.46 mg/kg in the 64 tested red varieties and from 1.36 to 30.21 mg/kg, with a mean of 10.83 mg/kg in the 36 tested white varieties. However, in our study, in the 2014 winter season berries of 'Victoria' grapes, the total flavonol concentration could be up to 100 mg/kg FW, which might owe to the high expression of VviFLS4 in the berries.

Effects of Growing Season on Berries Flavanol Composition
The influence of growing seasons on flavanols seemed to be variety dependent, which was different from the consistent influence found in the anthocyanins and flavonols accumulation in our study. Similarly, in a previous study, Zhu et al. [6] reported that the winter season berries of 'Khoyo' and 'Muscat Hamburg' had higher flavanol concentrations than the summer season berries at harvest, while no significant difference was found in the summer and the winter season berries of the 'NW196 grape. However, Xu et al. [7] showed that flavanols in the skin of the winter season berries were higher than those of the summer season berries for all of the cultivars they investigated. Although flavanols shared the same common upstream steps with flavonols and anthocyanins, the high temperature seemed to have a minor influence on the flavanol concentration in our study. Until now, the effect of temperature on the flavan-3-ol biosynthesis and PA accumulation was not well understood [45]. Some studies reported that no effect on the skin PAs when the high temperature treatments were applied at the whole-vine level [41,47], while other studies also reported a decrease in the skin PAs was correlated with higher temperature [48,49]. The formation of flavanols in skins mostly occurred in the early stage of the berry development, starting from the fruit-set with maximum levels observed before veraison [50]. Under the double cropping system, the summer season had less heat pressure and the winter season had a warm condition in the former stage of the berry development, which might cause little variation in the flavanol concentration between the summer and the winter season berries. VviLAR1 (VIT_201s0011g02960) and VviLAR2 (VIT_217s0000g04150) were upregulated in the 'Muscat Hamburg' winter season berries, which led to a higher flavanol concentration than the summer season berries in 2014. However, in the 'Victoria' grapes, VviLAR2 was downregulated in the winter season berries, which was in agreement with the higher flavanol concentration found in the 2014 summer season berries than those in the winter season berries in the previous analysis.

Conclusions
In this study, the cool weather conditions in the winter season made the grape berries accumulate more phenolic compounds than those in the summer season under the double cropping system. Although this cropping system could avoid intense rainfall and heatwave as much as possible, the high summer season temperature still inhibited the flavonoid metabolism to a great extent. For anthocyanins, the non-acylated and non-methylated anthocyanins were easily degraded in response to the high temperature, thus leading to the higher acylated and methylated anthocyanin proportions in the summer season berries. Flavanols were more stable to the climate changes than anthocyanins and flavonols, because the flavanol concentration did not show a significant decrease in the summer season berries. Most of the genes and transcription factors related to the flavonoid biosynthesis were upregulated in the winter season berries, which was consistent with the results found in the metabolites. The variation in VviLARs expression between 'Muscat Hamburg' and 'Victoria' grapes might be the reason why flavanols showed different trends in response to the climate variation among seasons. The extreme weather conditions in the summer season provide possible insights into how global warming or climate changes would impact viticulture in the future.

Data Availability Statement:
The data showed in this study are contained within the article.